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ABSTRACT 

We have made differential radial velocity measurements of the star a Cen A using two spectrographs, UVES 
and UCLES, both with iodine absorption cells for wavelength referencing. Stellar oscillations are clearly visible 
in the time series. After removing jumps and slow trends in the data, we show that the precision of the velocity 
measurements per minute of observing time is 0.42 ms~' for UVES and LOms"' for UCLES, while the noise 
level in the Fourier spectrum of the combined data is 1 .9 cm s"' . As such, these observations represent the most 
precise velocities ever measured on any star apart from the Sun. 

Subject headings: stars: individual {a Cen A) — stars: oscillations — techniques: radial velocities 



1. INTRODUCTION 

The search for extra-solar planets has driven tremendous im- 
provements in high-precision measurements of stellar differ- 

■ ential radial velocities (e.g., Marcy & Butler 2000). In par- 
, allel, the same techniques have been used with great success 

to measure stellar oscillations. Recent achievements, reviewed 
by Bouchy & Carrier (2003) and Bedding & Kjeldsen (2003), 
. include observations of oscillations in Procyon with ELODIE 

■ (Martic et al. 1999, see also Barban et al. 1999); (3 Hyi with 
; UCLES (Bedding et al. 2001) and CORALIE (Camer et al. 

2001); and a Cen A and B with CORALIE (Bouchy & Carrier 
2001, 2002; Carrier & Bourban 2003). Here we report obser- 
, vations of a Cen A made with UVES and UCLES which repre- 

■ sent the most precise differential radial velocities ever measured 
' on any star apart from the Sun. 

2. OBSERVATIONS 

We observed a Cen A in May 2001 from two sites. At the 
European Southern Observatory in Chile we used UVES (UV- 
Visual Echelle Spectrograph) at the 8.2-m Unit Telescope 2 of 
the Very Large Telescope (VLT)''. At Siding Spring Observa- 
tory in Australia we used UCLES (University College London 
Echelle Spectrograph) at the 3.9-m Anglo- Australian Telescope 
(AAT). In both cases, an iodine cell was used to provide a stable 
wavelength reference (Butler et al. 1996). 

At the VLT we obtained 3013 spectra of a Cen A, with typ- 
ical exposure times of 1-3 s and a median cadence of one ex- 
posure every 26 s. At the AAT we obtained 5169 spectra of 
a Cen A, with typical exposure times of 6 s and a median ca- 
dence of one exposure every 20 s. The resulting velocities are 
shown in Fig. 1, and the effects of bad weather can be seen (we 
were allocated four nights with the VLT and six with the AAT). 

The UCLES velocities show upward trends during most 
nights, which we believe to be related to the slow movement 
of the CCD dewar as liquid nitrogen boiled off. Nights three 
and (particularly) five also show a jump that coincides with the 
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refilling of the CCD dewar in the middle of the night. These 
UCLES observations differed from our earlier run on (3 Hyi 
(Bedding et al. 2001), and also from standard UCLES planet- 
search observations, in that the CCD was rotated by 90 degrees 
to speed up readout time. This had the side-effect of causing 
the CCD readout to be in the same direction as the dispersion, 
and also of making this direction vertical (so that flexure in the 
dewar due to changes in its weight would have shifted the spec- 
trum in the dispersion direction). It seems there was an effect 
on the velocities that was not completely corrected by the io- 
dine reference method, which suggests that either the PSF de- 
scription or the spectrum extraction from the CCD images was 
inadequate. We note that improvements to the PSF description 
are an active area of work to enhance immunity to spectrometer 
changes. The UVES data, meanwhile, show slow trends and 
two smaller jumps which are presumably also instrumental - 
at least one of the jumps can be identified with a correction to 
the position of the star on the slit. While these jumps and slow 
trends would seriously compromise a planet search, they fortu- 
nately have neglible effect on our measurement of oscillations. 

To remove the slow trends, we have subtracted a smoothed 
version of the data, and this was done separately on either side 
of the jumps so that they, too, were removed. The detrended 
time series are shown in Fig. 1. We have verified by calculating 
power spectra that this process of high-pass filtering effectively 
removes all power below about 0.2 mHz. It also removes any 
power at higher frequencies that arises from the jumps, which 
would otherwise degrade the oscillation spectrum. 

3. ANALYSIS AND DISCUSSION 

The time series of velocity measurements clearly show oscil- 
lations and the effects of beating between modes (Fig. 2). The 
data presented here have unprecedented precision and we are 
interested in obtaining the lowest possible noise in the power 
spectrum so as to measure as many modes as possible. We also 
wish to estimate the actual precision of the Doppler measure- 
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ments. 

The algorithm used to extract the Doppler velocity from each 
spectrum also provides us with an estimate cr; of the uncertainty 
in this velocity measurement. These are derived from the scat- 
ter of velocities measured from many small 1 .7 A) segments 
of the echelle spectrum. In the past, we have used these values 
to generate weights (w, = 1 / af) for the Fourier analysis (Bed- 
ding et al. 2001; Kjeldsen et al. 2003). In the present analysis, 
we take the opportunity to verify that these cj, do indeed reflect 
the actual noise properties of the velocity measurements. 

Our first step in this process was to measure the noise in the 
power spectrum at high frequencies, well beyond the stellar sig- 
nal. It soon became clear that this procedure implied a mea- 
surement precision significantly worse than is indicated by the 
point-to-point scatter in the time series itself. The implication is 
that some fraction of the velocity measurements are 'bad,' con- 
tributing a disproportionate amount of power in Fourier space. 

Since the oscillation signal is the dominant cause of varia- 
tions in the velocity series, we need to remove this signal in or- 
der to estimate the noise and to locate the bad points. We chose 
to remove the signal iteratively by finding the strongest peak 
in the power spectrum and subtracting the corresponding sinu- 
soid from the time series. This procedure was carried out for 
the strongest peaks in the oscillation spectrum in the frequency 
range 0-3.5 mHz, until spectral leakage into high frequencies 
from the remaining power was negligible. We were left with 
a time series of residual velocities, r,, which reflects the noise 
properties of the measurements. 

The next step was to analyse the residuals for evidence of 
bad points, which we would recognize as those values deviat- 
ing from zero by more than would be expected from their uncer- 
tainties. In other words, we examined the ratio rj/dj, which we 
expect to be Gaussian distributed, so that outliers correspond 
to suspect data points. We found that the best way to inves- 
tigate this was via the cumulative histogram of |?",7cr, |, which 
is shown for both the UVES and UCLES as the points in the 
upper panels of Fig. 3. The solid curves in these figures show 
the cumulative histograms for the best- fitting Gaussian distribu- 
tions. We indeed see a significant excess of outliers for \ri/ai\ 
beyond about 2, in both data sets. The lower panels show the 
ratio between the points and the curve, which is the fraction / 
of data points that could be considered as good. This fraction is 
essentially unity out to cr, ] ~ 2, and then falls off, indicating 
that about half the data points with | > 3 are bad. 

At this point, we could simply make the decision to ignore 
all data points with |r,/(T, | above a certain value, such as 3.0, 
on the basis that many of them would be bad points that would 
increase the noise in the oscillation spectrum. We instead chose 
a more elegant approach, which gave essentially the same re- 
sults, in which we used the information in Fig. 3 to adjust 
the weights: those points with large values of |?",7cr,| were de- 
creased in weight, with every w,- multiplied by the factor /. 
Given that the weights are calculated as w, = 1 /af, the adjust- 
ment was achieved by dividing each measurement error (tr,) by 
the square root of /. 

With these adjustments to the measurement uncertainties, 
which effectively down- weight the bad data points, we now ex- 
pect the noise floor at high frequencies in the amplitude spec- 
trum of the residuals (r,) to be substantially reduced. We mea- 
sured the average noise (Jamp in the range 7.5-15 mHz to be 
2.11cms-i for UVES and 4.37 cm s"! for UCLES. The cor- 
responding values before down-weighting the bad data points 



were 2.33 cm s"' for UVES and 4.99 cm s"' for UCLES. We 
can therefore see that adjusting the weights has lowered the 
noise by about 10%. 

The final stage in this processing involved checking the cali- 
bration of the uncertainties. By this, we mean that the estimates 
a, should be consistent with the noise level determined from the 
amplitude spectrum. On the one hand, the mean variance of the 
data can be calculated as a weighted mean of the cr;, as follows: 

N / N 

j2<^fwJY,^i, (1) 

1=1 / i=i 

where the weights are given by w, = 1 / af (which means the nu- 
merator is simply equal to A^). On the other hand, the variance 
deduced from the noise level cramp in the ampHtude spectrum is 
(Appendix A. 1 of Kjeldsen & Bedding 1995): 

crLpN/TT. (2) 
We require these to be equal, which gives the condition 

N 

'^lap^'^f = '^- (3) 

!=1 

Using the values of (Tamp given above, we concluded that Eq. (3) 
would be satisfied for each data set provided the uncertainties 
fT, were multiplied by 0.78 for UVES and 0.87 for UCLES. It 
is these calibrated cr, that are shown in Figs. 2 and 4, and they 
represent our best estimate of the high-frequency precision of 
the data. Also shown in Fig. 4 as solid lines are the residuals 
r, after smoothing with a running box-car mean. We can see 
that there is generally very good agreement between these two 
independent measures of the velocity precision, giving us confi- 
dence that we have correctly estimated cr, both in relative terms 
and in the absolute calibration. 

We have also investigated the dependence of the velocity pre- 
cision on the photon flux in the stellar spectrum. Figure 5 shows 
cr, versus the signal-to-noise ratio (SNR), where the latter is the 
square root of the median number of photons per pixel in the 
iodine region of the spectrum. The points for both UVES and 
UCLES agree well with a slope of -1 in the logarithm, as ex- 
pected for Poisson statistics. We also note that the offset be- 
tween the two sets of points arises because of the difference in 
dispersion of the spectrographs. UVES has higher dispersion, 
by a factor of 1.55, which means that not only is velocity preci- 
sion expected to be better by this factor, but also that the SNR 
per unit wavelength is lower (by the square root of this factor). 
Combining these effects, we would expect the two distributions 
to be in the ratio 1.55''^ = 1.93, which is exactly what we find 
in Fig. 5. In other words, the lower precision in the UVES data 
at a given SNR per pixel results from the higher dispersion but, 
allowing for this, the intrinsic precision is the same for the two 
systems. In addition, the CCD on the UVES system was able 
to record more photons per exposure than UCLES (the median 
SNR is considerably higher). 

The result of the processing described above was a time se- 
ries for both UVES and UCLES, each of which consisted of the 
time stamps, the measured velocities (after correction for jumps 
and drifts) and the uncertainty estimates (after the adjustments 
described). These two time series could then be merged in order 
to produce the oscillation power spectrum, and this is shown in 
Fig. 6. The average noise in the ampUtude spectrum in the range 
7.5-15mHz is 2.03 cms"'. As discussed at the beginning of 
this section, some of this power comes from spectral leakage of 
the osciUation signal itself. Therefore, a more accurate measure 
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of the noise is obtained from the power spectrum of the residu- 
als, r,-, when using the final weights, l/crf. The result gives an 
average noise in the range 7.5-15 mHz of 1 .9 1 cm s~' . For com- 
parison, the noise level reported from CORALIE observations 
of a Cen A was 4.3 cms"' (Bouchy & Carrier 2002), while 
observations of a Cen B gave 3.75 cms"' (Carrier & Bourban 
2003). 

We can also calculate the precision per minute of observing 
time, which is shown in Table 1 . For comparison, we include 
velocity precision from other oscillation measurements. The 
list is not meant to be exhaustive, but most of the instruments 
that have been used in recent years are represented. The preci- 
sion depends, of course, on several factors such as the telescope 
aperture, target brightness, observing duty cycle, spectrograph 
stability and method of wavelength calibration. It is clear that 
the observations reported here, particularly those with UVES, 
are significantly more precise than any previous measurements 
of stars other than the Sun. Of course, we refer here to the pre- 
cision at frequencies above ~ 0.8 mHz, which is the regime of 
interest for oscillations in solar-type stars. 

The referee has questioned whether adjusting the weights us- 
ing the method described above might have affected the accu- 
racy with which the oscillation frequencies can be measured. 
To test this, we have extracted the ten highest peaks from both 
the power spectrum in Fig. 6 and from the power spectrum ob- 
tained without adjusting the weights. The frequencies of all ten 
peaks agreed very well, with a mean difference of 0. 11 /uHz and 
a maximum difference of 0.4 /xHz. The latter is a factor of ten 



smaller than the FWHM of the spectral window and probably 
well below the natural linewidth of the modes. Therefore, there 
is no reason to think that the reduction in the noise level ob- 
tained by adjusting weights has come at the expense of reduced 
accuracy in the measured frequencies. 

4. CONCLUSION 

We have analysed differential radial velocity measurements 
of a Cen A made with UVES at the VLT and UCLES at the 
AAT. Stellar oscillations are clearly visible in the time series. 
Slow drifts and sudden jumps of a few metres per second, pre- 
sumably instrumental, were removed from each time series us- 
ing a high-pass filter. We then used the measurement uncer- 
tainties as weights in calculating the power spectrum, but we 
found it necessary to modify some of the weights to account 
for a small fraction of bad data points. In the end, we reached 
a noise floor of 1.9 cms"' in the amplitude spectrum and in a 
future paper we will present a full analysis of the oscillation 
frequencies and a comparison with stellar models. 

We thank Bill Chaplin for providing a times series that al- 
lowed us to estimate the precision of BiSON observations. 
This work was supported financially by the Australian Research 
Council, by the Danish Natural Science Research Council and 
by the Danish National Research Foundation through its estab- 
lishment of the Theoretical Astrophysics Center. We further 
acknowledge support by NSF grant AST-9988087 (RPB), and 
by SUN Microsystems. 
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Fig. 1 . — Time series of velocity measurements of a Cen A, both before and after removal of jumps and slow trends. The times of the jumps in UVES are at 
6.655 and 6.760, and those for UCLES are at 6. 122 and 8. 1 10. 
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Fig. 2. — A four-hour segment of the detrended UVES velocity time series, showing l-cr error bars. 
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Fig. 3. — Upper: cumulative histograms of for UVES (left) and UCLES) riglit. The points show the observations and the soHd curve shows the result 

expected for Gaussian-distributed noise. Lower: the ratio between observed and expected histograms, indicating the fraction of "good" data points. 
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Fig. 4. — Points show measurement errors in each time series, while the (red) curves show the smoothed residuals (see text) 



Table 1 

Noise levels from observations of solar-like oscillations 



Star 


Spectrograph 


Precision in 
1 min (m s"' ) 


Reference 


Sun 


BiSON 


0.2 


data supplied by W. Chaplin 


a Cen A 


UVES 


0.42 


this paper 


Sun 


GOLF 


0.6 


data from the GOLF Web site 


a Cen A 


UCLES 


1.0 


this paper 


aCenB 


CORALIE 


1.7 


Carrier & Bourban (2003) 


a Cen A 


CORALIE 


1.7 


Bouchy & Carrier (2002) 


Procyon 


ELODIE 


2.5 


Martic et al. (1999) 


Procyon 


CORALIE 


2.7 


Carrier et al. (2002) 


/3Hyi 


UCLES 


3.0 


Bedding et al. (2001) 


/3Hyi 


CORALIE 


4.2 


Carrier et al. (2001) 


Procyon 


AFOE 


4.2 


Brown (2000) 


Procyon 


FOE 


4.7 


Brown et al. (1991) 


Procyon 


HIDES 


5.1 


Kambe et al. (2003) 
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Fig. 5. — Measurement errors as a function of SNR. The two straight lines both have slopes of -1, and are displaced from each other by a factor 1.93 (see text). 
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Fig. 6. — Power spectrum of the combined velocity time series. The inset shows the spectral window, with the frequency scale expanded by a factor of ten. 



